NiO - from first principles 
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We have calculated from first-principles the octupolar interactions of the Ni 2+ ion in NiO, which 
gives the theoretical basis for the ionic description of properties of NiO with fully localized strongly- 
correlated eight d electrons. A failure of the up-now first principles ionic calculations for NiO was 
largely due to too small values taken for the octupolar moment of the transition-metal atom, largely 
generated by too small value for (r$). Our many-electron crystal-field based approach enables 
successful calculations of the electronic structure and magnetic properties both in the paramagnetic 
and in magnetically-ordered state as well as zero-temperature properties and thermodynamics. 
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During last 10 years 0, 0, IS 0] we have described 
the magnetism and low-energy electronic structure of 
NiO within the localized atomistic paradigm with only 
three parameters: the octupolar crystal-field parameter 
B± =+21 K (10Dq=1.086 eV), the spin-orbit coupling 
A 5 ^, =-480 K and a small trigonal distortion B\ =+50 
K [5j. For the magnetic state the fourth parameter, the 
molecular-field coefficient n, has been introduced - n — 
—200 T/hb yields T/v of 525 K in agreement with the 
experimental observation. The fundamentally important 
value of B4 has been fitted to the excitation of 1.06-1.13 
eV seen in the optical absorption 0] . Here we would like 
to calculate £4 from first principles - it buckles the inter- 
nal consistency of our understanding of NiO, explaining 
its insulating ground state, magnetism both in the AF 
and paramagnetic state. 

Although the ionic crystal-field based approach seems 
to be the most natural approach to transition-atom ox- 
ides and is known for 75 years starting with pioneering 
works of Bethe and Van Vleck the presently in fash- 
ion are theoretical approaches pointing out the itinerant 
(band) character of d electrons 0, S EJ EH 13 US 13 
This paper has been motivated by a series of ab initio and 
first-principles papers devoted to explanation of proper- 
ties of NiO, which recently have been published within 
the ba nd p icture and the continuous energy spectrum 
EEinfEHIIl- We note that with time the orig- 
inal LDA approach is significantly modified by incorpo- 
ration of different local potentials, i.e. goes towards the 
advocating by us many-electron crystal-field approach. 

From historical point of view, the band description has 
been introduced as an opposite view to the ionic view of 
Van Vleck, which underlies the discrete electronic struc- 
ture. After works of Van Vleck a milestone achievement 
within the localized paradigm were the Tanabe-Sugano 
diagrams known from year of 1954 0,H2 an( A shown in 
Fig. 1 for the Ni 2+ ion. A quite successful description in 




FIG. 1: A Tanabe-Sugano diagram for the Ni ion (3d con 



figuration) showing the effect of the octahedral crystal field 
on the electronic terms of the free Ni 2+ -ion. The electronic 
structure of cubic subterms, corresponding to Dq/B =0.85, 
relevant to NiO is marked by an arrow. Zero-field term values 
are presently well known in the NIST database. 

1984 by Fujimori and Minami [2^] of the photoemission 
data by a local configuration-interaction cluster model 
was very soon overflowed by the band-based view point- 
ing out the dominant role played by the covalency and 
strong hybridization of Ni 3c? and O 2p states. 

In the present situation of the dominant band view 
in the NiO problem, in general in 3d oxides, we feel 
necessary to complete our understanding of properties 
of the magnetism and electronic structure of NiO in 
the localized paradigm. Our understanding is based 
on well-known physical concepts like the crystal-field 
(CEF) |24j, spin-orbit (s-o) coupling, local distortions 
and other terms known from the ionic language. We 
present this view being aware that the ionic picture and 
crystal-field considerations are at present treated as the 
"old-fashioned" and contemptuous physics in times of 
wide spreading omnipotent band theories of different ver- 
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FIG. 2: The calculated temperature dependence of the 3 F 
term (Ni 2+ , 3d 8 ) in NiO in the presence of the octahedral crys- 
tal field and spin-orbit coupling in the magnetically-ordered 
state below Tjv of 525 K. In the paramagnetic state the elec- 
tronic structure is temperature independent unless we do con- 
sider, for instance, a changing of the CEF parameter due to 
the thermal lattice expansion. The states of the 3 A2 g subterm 
are characterized by magnetic moment of 0, +2.26 and -2.26 
fiB- In the magnetic state there occurs the spin polarization 
of this three-fold degenerated 3 A2 g subterm. A small split- 
ting, 1.0 meV, of the 3 ' Ai g subterm in the paramagnetic state 
due to the presence of the trigonal distortion is not visible in 
this energy scale. 

sions LDA, LSDA, LDA+U % EH, LDA+GGA [l^, 
LDA+DMFT |20| and many, many others. We gave a 
name of QUASST for our approach to a solid containing 
transition-metal atoms from Quantum Atomistic Solid 
State Theory pointing out that the physically adequate 
description of properties of a 3d/4f/5f solid the best is to 
start from analysis of the electronic structure of consti- 
tuting atoms. In a solid atoms occur rather as ions by 
which we understand a well-defined charge state of the 
given atom. 

Schematic steps of the QUASST approach to NiO can 
be written as: 

1. we accept the atomistic structure of matter, what 
means that atoms preserve much of their atomic prop- 
erties becoming the full part of a solid. In NiO during 
the formation of the compound there occurs an electron 
transfer of two electrons from Ni to O with the simul- 
taneous formation of the rocksalt NaCl lattice with the 
ionic charge distribution Ni 2+ 2_ ; 

2. all magnetism of NiO is related to the Ni 2+ ion with 
eight electrons outside close configuration 18 Ar because 
2 ~ ions have closed shells ( 10 Ne); 

3. eight electrons of the Ni 2+ ion form strongly- 
correlated atomic-like system, 3d 8 ; its electronic struc- 
ture is known from the atomic physics (NIST 25]: l D - 
1.74 eV, 3 P - 2.10, l G - 2.86; l S - 6.51 eV with the 3 F 
Hund rule ground term); 

4. from the translational symmetry of the rocksalt 



structure of NiO we learn that a) all Ni sites are equiv- 
alent, b) the local symmetry of the Ni cation is octahe- 
dral, c) the Ni ion is surrounded by six nearest neighbor 
oxygens forming almost perfect octahedron, d) with the 
lowering temperature there appears small off-octahedral 
trigonal distortion; 

5. influence of the octahedral crystal field on the free- 
ion electronic terms as is shown in the Tanabe-Sugano 
diagram, Fig. 1; 

6. influence of the intra-atomic spin-orbit coupling on 
the local electronic structure; for the free Ni 2+ ion A = 
-42 meV; NIST @ gives the multiplet splitting 3 F - 0.0 
( 3 F 4 ), 168.7 meV ( 3 F 3 ) and 281.4 meV ( 3 F 2 ); 

7. influence of the off-octahedral rhombohedral distor- 
tion on the local electronic structure; 

8. magnetic interactions lead below 525 K to the 
magnetically-ordered state; 

9. having determined the local electronic structure 
with the eigenfunctions we have the zero-temperature 
properties as well as the free energy -F(T); 

10. the zero-temperature properties are related to the 
properties of the local ground state; the free energy F(T) 
enables calculations of the whole thermodynamics. See, 
for instance, our description of FeBr 2 0, an exemplary 
of the magnetic 3d compound. 

Ad. 1. An assumption about the atomistic construc- 
tion of matter seems to be obvious but we would like to 
point out that band-structure calculations for 3d oxides 
disintegrate 3d atoms completely starting from consid- 
eration of all 3d electrons as independent electrons in 
the octahedral crystal field. According to us these one- 
electron calculations have to reproduce, before calcula- 
tions of properties of a solid, the electronic structure 
of the given ion, i.e. its term structure. For NiO, the 
3 A 2g , 1 E and 1 A\ subterm structure of the 6-fold degen- 
erated t%g e2 g configuration. In our approach we simply ac- 
cept the term electronic structure known from the atomic 
physics (see Ad. 3). We point out that the intraatomic 
term structure results from strong electron correlations 
among d electrons. 

Ad. 3. From the atomic physics we know that for the 
3d 8 configuration of the Ni 2+ ion 45 states are grouped 
in 5 terms: 3 F (21 states), 3 P (9), 1 G (9), 1 D (5), 1 S 
(1). The first excited term is almost 2 eV above being 
inactive for low- and room-temperature properties. 

Ad. 4. For more complex structures there will be a 
few 3d sites, each of them having own electronic struc- 
ture. Analysis of low-temperature structure is important 
to decide if all Ni sites are equivalent and all ions equally 
contribute to macroscopic properties. 

Ad. 5. Influence of the octahedral crystal field on 
the free-ion electronic terms has been calculated by Tan- 
abe and Sugano in a year of 1954 already |2jJ,|22j. The 
splitting of electronic terms to octahedral subterms in a 
function of the strength of the octahedral CEF param- 
eter Dq/B (B - intra-ionic Racah parameter determines 
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the energy scale) is known as Tanabe-Sugano diagrams, 
Fig. 1. A problem was and still is with i) the accep- 
tance of its validity for a solid compound and with ii) 
the evaluation of the value of Dq/B on this diagram for 
a given compound. Numerous qualitative indications, 
starting already at fifties of the XX century, have not 
been conclusive. In the crystal-field theory parameter 
lODq (=600-i?4) in the simplest form is the multiplica- 
tion of the octupolar charge moment of the lattice A4 and 
of the octupolar charge moment of the involved cation 
caused by anisotropic charge distribution of the own in- 
complete 3d shell. Thus Dq (B4) can be calculated from 
first principles provided the octupolar charge moment, /3 
(r 4 .), of the involved ion is known. (3 is the fourth-order 
Stevens coefficient. 

The octahedral crystal field coefficient A4 that is the 
octupolar charge moment of all surrounding charges at 
the Ni site can be calculated from the point-charge model 
which is first-principles elementary calculations. Taking 
the charge of oxygen Z as -2|e| and the nearest octahedron 
with the cation-oxygen distance d of 208.5 pm in NiO we 
obtain by formulae A 4 =7-Ze 2 /16-d 5 [24| a value for A4 of 
+290 Ka^ 4 , a s is the Bohr radius. Taking for the Ni 2+ 
ion f3 — +2/315 and a widely-spread Hartree-Fock result 
by Freeman- Watson (F-W) (r 4 ) =3.003 a.% a value 
of B4 of +5.5 K is obtained only. This value is 3.8 times 
smaller than the recent experimental determination of 
B4 of +21 K In our papers for 3d oxides (LaMnC>3, 
FeBr 2 , CoO, ...) we pointed out that despite this differ- 
ence the most important is that these ionic ab initio cal- 
culations give the proper sign of the B4 parameter - it is 
important because it determines the ground state of the 
cation and the experimentally derived strength of crystal- 
field interactions turns out to be much weaker than it was 
thought so far in literature for justification of the strong 
crystal-field approach in which the one-electron approach 
becomes more physically adequate. Here we would like 
to explain this theoretical 3.8-times smaller value - we 
attribute it to a large underestimation of (r 4 .) =3.003 a B 
obtained in F-W Hartree-Fock calculations of in a year 
of 1965 [H. 

At first, we have calculated the contribution to the oc- 
tupolar lattice potential from next neighbors - it strongly 
decreases due to the increase of the Ni-0 distance. The 
next-nearest contribution, due to the Ni 2+ ions, amounts 
to 9% of the first-octahedron contribution but it adds 
making theoretical evaluation of A4 as +315 Ka^ 4 . Thus, 
we claim that the proper value of (r 4 .) for the Ni 2+ ion 
is 10.50 a B - it is 3.5 times larger value than the F- 
W Hartree-Fock result from 1965. We claim that this 
widely-spread early Hartree-Fock result is erroneous. In- 
deed, the presently derived values are significantly larger. 
Korotin et al. [27j have used a value for of 1.26 A yield- 
ing approximately (r 4 .) = 32.2 afj. Recently Solovyev 
[28j, page 5, has calculated for the Ti 3+ ion (r 2 ,) of 2.27 
A 2 (=8.11 d? B ) : which approximately yields (r 4 .) even 
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FIG. 3: Comparison of the calculated electronic structure, 
shown in Fig. 2, with that experimentally observed by means 
of optical absorption in Ref. ||j . 

as large as 65 &%. Thus we think that a value of 10.5 
a|j needed for the purely ionic electrostatic-origin of the 
crystal-field splitting in NiO is fully reasonable. 

Ad. 6. It is fundamentally important to take the spin- 
orbit coupling into account despite its weakness. Its con- 
sequence is that the physically adequate space becomes 
the spin-orbital space not the orbital space separated 
from the spin space. It generates, for instance, physi- 
cally adequate electronic structure of d electrons and the 
orbital moment. 

Ad. 7. The rhombohedral (trigonal) distortion is small 
and causes a slight splitting of the lowest quasi-triplet 
in NiO with D of order of a few meV (for £ 2 value of 
+50 K shown in Fig. 2 D amounts to 1.0 meV only). 
But this small distortion determines in our calculations 
the direction of the Ni magnetic moment, in NiO the 
(112) direction, i.e. perpendicularly to the trigonal axis. 
Such moment direction causes further symmetry lowering 
due to magnetostriction. Thus, we can say that we fully 
understand, and we calculate, a delicate interplay of the 
magnetism of NiO and distortions of its crystal structure. 

Ad. 8. The formation of the magnetic state we 
have described in numerous compounds - let mention ex- 
emplary 4f/3d/5f compounds ErNi 5 , FeBr 2 , NiO, CoO, 
UPd2Al 3 and UGa2, results of which have been published 
starting from 1992. In all these cases the magnetic energy 
is much smaller than the overall CEF splitting. They are 
both ionic (FeBr 2 @, NiO, CoO 0) and intermetallic 
(ErNi 5 , UPd 2 Al 3 , UGa 2 ) compounds. 

We would like to note that all of the used by us pa- 
rameters (dominant octahedral CEF parameter B4, the 
spin-orbit coupling A S _ Q , lattice distortions) have clear 
physical meaning. The most important assumption is the 
existence of very strong correlations among 3c? electrons 
preserving the atomistic ionic integrity of the Ni 2+ ion 
also in the solid NiO. We stress the good reproduction 
of experimental results like temperature dependence of 
the heat capacity 0, the value of the magnetic moment 
and its (112) direction. Our calculations yield a total 
magnetic moment of 2.53 \ib in very good agreement 
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with recent experimental evaluation |3lJ, |32j. The cal- 
culated moment contains almost 20% contribution from 
the orbital moment. Our purely ionic crystal-field ap- 
proach, despite its simplicity, is in agreement with recent 
experimental findings of Jauch and Reehuis |33|. who 
have revealed by the high-accuracy single crystal 7-ray 
diffraction experiment the electron-density distribution 
in NiO indicating i) the absence of the covalency and ii) 
the Ni-0 interaction to be purely ionic. For photoemis- 
sion results we largely accept the Fujimori-Minami inter- 
pretation [23j as their cluster result about the multiplet 
structure, shown in their Fig. 7, is quite similar to our 
purely ionic multiplet structure, shown in Figs 1 and 3. 

In conclusions, we have calculated from first- 
principles the octupolar interactions of the Ni 2+ ion 
in NiO, which gives the theoretical basis for the ionic 
description of properties of NiO with fully localized 
strongly-correlated eight d electrons. The crystal-field in- 
teractions are relatively strong in NiO but not so strong 
to destroy the ionic integrity of the 3d electrons. The in- 
sulating gap in NiO of 4.3 eV is of the same origin as a gap 
in classical insulators, MgO or CaO. A failure of the up- 
now first principles ionic calculations for NiO were largely 
due to too small values taken for the octupolar moment 
of the Ni 2+ ion, largely generated by too small value for 
(rj). Our many-electron crystal- field based approach en- 
ables successful calculations of the electronic structure 
and magnetic properties both in the paramagnetic and 
in magnetically-ordered state as well as zero-temperature 
properties and thermodynamics. Thus we claim that 
this many-electron CEF picture is the physically ade- 
quate starting point for description of magnetism and 
low-energy electronic structure of 3d oxides. Paraphras- 
ing conclusion of Ref. [2(| we can say, that in spite of the 
limitations of the implementation of magnetism and spin- 
polarization we showed that the ionic CEF-based ap- 
proach which combines first-principles, material-specific 
information with very strong electron correlation effects 
is able to deal with late transition-metal monoxides like 
NiO. This old-known compound turns out to exhibit sub- 
picosecond switchin g ch aracteristics of potential applica- 
tions in spintronics |34l|. 
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